
Mad# tvaflable under NASA sponsorship 
in thejnteresi of early and wide dis- 
semination of Earth Resources Survey 
Program information and without liability 
for any use made thereof.” 


Bimonthly Report: 11/4/73 - 1/4/74 


ERTS Proposal No, 108 
Remote Sensing of Ocean Currents 


E 7* 4 ~ 1 0* ^ 1 9* 

fOp 


George A. Maul, Principal Investigator, GSFC ID-C0315 


43- 

O 

i— 

V3 

r- 

<d Ch 

f • 

H *■" 

& 

o 


o 

55 

& o 

j; 

i 





V 

| 

no 

O 

1 

j 

U 


00 

■ 

o 

f 

n 


r- 

u 

1 7) 

to 

Ss r- 

u o 



w • a 
U > aJ 

CL 

o o 

SS U fN 

Pm •r“t rn 

o a 

<d 

o % a) 


S3 -P u 


H M o 

o 

to o 

•H 

SH +• 

w a> cd 

(d 

wp 5 a 

M 

!• o 

-H 


(/I 

' Eh H -P 


q 43 «d 

G 

53 -M S3 


w a ^ 

s 

M Q 


1 s 

< 

! -H r* 

o 

i;ct> v~ 

•ri 

!i- 

U 

m <o • 

<D 

o eh a 

!r- sa <d 

O- 

;l W ^ 

w in 

j r* C* 

o ^ 

r- « =r 

e • 

i« o 


U 1 

^3 tO 


OBJECT 

The object of this investigation is to locate ocean current boundaries by 
sensing the color change associated with the cyclonic edge of the zone of 
maximum horizontal velocity shear. The test site is the eastern Gulf of 
Mexico where the strongly baroclinic flow from the Yucatan Straits forms 
into the Loop Current. The research will attempt to use ERTS data in the 
investigation of ocean color sensing from simultaneous observations by ship 
and satellite. 


FIELD DATA COLLECTION 


The time series was completed with the September 1S73 cruise. Only short 
trips were performed during the reporting period for the purpose of observing 
spectra in affiliation with other ocean color experiments (one with GSFC and 
one with SKYLAB). The 1/4 meter Ebert spectroradiometer was lost from the 
project and further measurements are severely curtailed and depend on bor- 
rowed equipment. An effort is being made to purchase a new instrument. 


WORK SUMMARY 


Routine data processing and end of the year leave consumed most of the project 
time for this period. A paper entitled "Relationships between ERTS radiances 
and Oceanic Fronts’', by George A. Maul and Howard R. Gordon was presented 
by Maul at the 10 December 1973 Symposium on significant results in Washington, 
D.C. The authors are frantically trying to ready a manuscript before the 
publication deadline. 


WORK PLANS 

Continued studies on the theory of upwelling spectra and the relationship to 
ERTS radiances will be carried out. A SKYLAB experiment to supplement the 
observations is planned for 8 January. The text of a paper on computer en- 
hancement of ERTS scenes for the oceans is to be started. 


ghategraphy may ba 
iitaf Dsta Center 
iifth Dakota Avinu§ 
mm Falls,. SD 





Reproduced from 
best available copy. 





RELATIONSHIPS BETWEEN ERTS 
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ABSTRACT 


A time series of the Loop Current in the 
Gulf of Mexico, covering an annual cycle of 
growth, spreading, and decay, has been ob- 
tained in synchronization with ERTS. Com- 
puter enhanced images, which are necessary 
to extract useful oceanic information, show 
that the current can be observed either by 
color or sea state effects associated with 
the cyclonic boundary. The color effect 
relates to the spectral variations in the 
optical properties of the water and its sus- 
pended particles, and is studied by radi- 
ative transfer theory. Significant oceanic 
parameters identified are: the probability 

of forward scattering, and the ratio of 
scattering to total attenuation. Several 
spectra of upwelling diffuse light are com- 
puted as a function of the concentration of 
particles and yellow substance. These 
calculations compare favorably with experi- 
mental measurements and show that the ratio 
of channels method gives ambiguous inter- 
pretative results. These results are used 
to discuss features in images where surface 
measurements were obtained and are extended 
to tentative explanation in others. 



DISCUSSION 


A copy of the paper presented at the symposium is attached and will serve as 
a synopsis of our current thoughts. One point about the symposium should be 
brought out: many ocean-oriented investigators do not seem to be aware of 

the physical principles involved in light from the sea, and it appears that 
some feedback mechanism should be initiated so that we all appreciate the 
problems better. Perhaps some sessions could be planned at a symposium in 
the future devoted to the mundane task. 


George A. Maul 
NOAA/AOML/PhOL 
15 Rickenbacker Causev/ay 
Miami, Florida 33149 



INTRODUCTION 


Remote sensing o'f the ocean in the visible region of the 
spectrum is being explored for a variety of reasons. Maximum 
insolation at 475 nm provides a natural energy source for 
passive sensors. A minimum in the attenuation coefficient 
for water may allow a passive measurement of some oceanic 
properties as a function of depth. Variations in the spectra 
of upwelling light at these wavelengths can be attributed to 
variations in pigment forming molecules such as in phyto- 
plankton and to variations in the concentration of scattering 
particles such as suspended sediments. Strongly baroclinic 
currents in the subtropics and tropics can be detected by 
change in the optical properties across their boundaries when 
infrared techniques fail due to isothermal surface conditions. 
Patterns of man’s activities in and on the ocean have near 
surface manifestations which change the nature of upwelling 
visible radiance. 

The social significance of oceanic observations from a 
vehicle such as ERTS has recently been summarized by Maul 
(1974a). Discussions in the present paper are limited to a 
study of the dominant feature of the Gulf of Mexico circula- 
tion, that portion of the Gulf Stream System called the Loop 
Current. This intense current transports vast amounts of 
thermal and kinetic energy into the basin through the Yucatan 
Straits. Temporal and spatial variability in the flow (figure 1) 
was studied as part of the ground truth time-series which 
followed the anti -cyclonic turning from Yucatan to the Florida 
Straits. Some dynamical results were reported by Maul (1973a) 
and will not be discussed here. As a purely descriptive 
oceanographic study however, these data have significance in 
hurricane intensification studies such as by Leipper and 
Volgenau (1970); the investigation of the Florida east coast 
red tide (Murphy, et at., ;i973) used these ship tracks and a 
supporting ERTS image to document their arguments; the en- 
vironmental impact of a recent accidental jettisoning of 
cyanide canisters following a ship collision in the Gulf of 
Mexico is being studied using these data (Corwin and Richardson, 
1974). Continued reconnaissance of this current pattern is 
significant in our understanding of, and reaction to, the 
transport, which affects fishing industries, marine trans- 
portation, and public health. 


ERTS OBSERVATIONS OF OCEANIC FRONTS 

The location of the cyclonic boundary of the Loop Current 
was obtained by tracking the 22 °C isotherm at 100 meters depth 
using expendable bathythermographs. The pathline of this 
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isotherm is of the order of 20 kilometers to the right 
(facing downstream) of the surface frontal zone (Hansen and 
Maul, 1970). This, coupled with ship observations of chloro- 
phyll- a, surface temperature, volume scattering function, and 
sea state, provides the baseline measurements to insure that 
the interpretations to follow are well founded. Figure 2 is a 
temperature-depth profile across the current along a suborbital 
track during a satellite transit; surface observations are at 
the top of the figure. This summarizes the general conditions 
of increasing temperature and salinity and decreasing chloro- 
phyll and scattering when • crossing the boundary into the 
current and serves to orient them to the indicator isotherm. 
These observations confirm (Maul and Hansen, 1972) that 
changes in all these properties occur simultaneously allowing 
an increase in confidence when a recognition decision is made 
using several variables. 

Diffuse reflectance from beneath the ocean, which is de- 
fined as the ratio of upwelling to downwelling irradiance, 
measured just above the surface, is rarely more than 0.05. 
Reflectance from the ocean’s surface, which is independent 
of diffuse reflectance, can be comparable or even substantially 
larger depending on sea. state, and this has been shown by Maul 
(1973b) to be a useful indicator of the current. Reflectance 
from clouds and agricultural scenes is sometimes an order of 
magnitude greater than from the ocean, even in the 500-600 nm 
wavelength (a.) region. In order for the NASA Data Processing 
Facility (NDPF) to produce an image for an average scene 
radiance, the ocean signal is compressed into the lowest few 
gray scales. This is clearly illustrated in figure 3 which is 
a scanline plot across the boundary of the Loop Current from 
the multispectral scanner. The large spikes in all four 
channels (MSS 4 upper) are clouds; there seems at first glance 
to be very little change in digital number (DN) , which is pro- 
portional to radiant intensity, as a function of the sample 
number. However, careful examination shows that the average 
value of the DN at samples greater than number 950 is slightly 
larger than those before this point. It will be seen that 
this marks the transition to higher radiances caused by in- 
creased sea state in the current. 

In order to graphically display this small change over a 
two-dimensional region, computer enhancement is necessary. 
Contrast stretching for the ocean scene is accomplished by 
first studying the frequency distribution of DNs in a training 
area on the image; the area is selected to be representative 
of the ocean away from land and is large enough to be statisti- 
cally significant. As seen from figure 3, such a histogram 
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would be strongly bimodal due to clouds. Selection of the upper 
(DN ) and lower (DN n ) cutoff can be set at + 2 standard deviations 
about the ocean mode; in this case 9<£N<13 was used. It should 
be noted that each image will have different cutoff values and 
that the final image is quite sensitive to the range of DN chosen 
(Charnel 1 et at., 1973). The general expression for computing 
the stretch variable (£ ) is; 


C = M 


DNu-DN 


DN u -DfM| 


n 


(for DNj<DN<DN u ) 


CD 


where M is the maximum value allowed by the output device, and 
n is an arbitrary integer exponent. Equation 1 produces a nega- 
tive of the input digital image; positive whole integers n, 
further stretch the low radiance values encountered m the ocean. 
The graphic result of using equation 1 on the data from which 
figure 3 was taken is given in figure 4. 

Figure 4 is a negative print of an area due north of the 
Yucatan Straits using MSS 5 data. Computer enhancement in this 
iamge uses only five gray scales of the 128 levels available, 
all values below DN-. are set to 127 and all above DN U are set to 
0 The boundary between the resident Gulf waters (left) and the 
current (right) is seen as a transition from light to dark tones 
respectively. Since the radiance levels m the ocean are so low. 
the oceanographer must resort to computer enhancement as this 

example shows. 

Figure 5 is another enhancement of the Loop Current boundary 
using MSS 5 where 7^DN<15 and n = 2. In this negative image of 
the western Florida Keys, water fTom Florida Bay extends into 
the Florida Straits and is entrained by the Loop Current . The 
current boundary in both figures 4 and 5 was delineated by sur- 
face vessel tracks during the day of the satellite transit . 

Notice that the current is darker in tone (higher m radiance) 
in figure 5. This is caused by the dominance of surface reflec- 
tion due to higher sea state in the current in figure ».j S !E 0 ^ 
pared to higher reflectance due to particles in the Florida Bay 

water in figure 5. 


SOME THEORETICAL CONSIDERATIONS 

At this point it is useful to consider in detail the proc- 
esses contributing to the radiance spectrum N(a) at the position 
of the satellite. Solar radiation incident at the top of the 
atmosphere is absorbed and scattered in the atmosphere. ome 
radiation is scattered back into space without striking the 
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ocean contributing a radiance N (a) at the satellite. The rest 
(which is not absorbed) will interact with the ocean. This 
interaction can yield upwelling radiance above the ocean m three 
ways: 1) specular reflection from the surface; 2) diffuse re- 

flection from foam (bubbles) on or just beneath the ocean surface, 
and 3) the diffuse reflection from water molecules and suspended 
particles in the water. Of these the first two phenomena are 
closely related in that they depend on the sea state, while the 
third source of radiance is essentially independent of sea state. 
The specular reflectance from the rough ocean siirface can be 
computed by the methods of Cox and Munk (1954) in terms of t e 
wind speed and incident radiance distribution. The diffuse re 
flectance from white caps can be approximately accounted for by 
assuming they are "white”, Lambertian, and have an albedo of 1 
so that they contribute uniform upwelling radiance just above 
the surface given by 

H o IF 


where f is the fraction of the scene covered by the white caps 
and H n (A) the irradiance incident on the sea surface. The third 
source of radiance, that from beneath the surface, is the mos ™ 
difficult to compute and will be discussed in detail below. We 
can write the radiance at the satellite as 


N(X) = N S (X) + yU)N ss (X) + a(X) N d (X) 


( 2 ) 


where N C ^(A) is the contribution at the surface due to reflec- 
tion from the surface and white caps, Nj is the diffuse radiance 
-just above the surface due to photons that have penetrated the 
ocean, a (A) and v(a) are atmospheric transmittance factors for 
Nj and N cc . « (a) and 7 (a) are_in general not equal since 


m 


Nj and N«. s . « (a) and 7(a) are m general not equal since 

the radiance distribution (variation with angle) of N ss and N, 
are different. It is possible for photons to be reflected frB 
the surface, backscatter from the atmosphere into the ocean, 
and scatter back into the atmosphere. Photons which do this we 
consider to be a part of Nj. It should be stressed at this 
point that Nj (A) in the above equation is the only source of 
radiance that contains information about conditions beneath the 
sea surface such as the concentration and composition of sus- 
pended particles and dissolved organic material. A thoroug 
understanding of the dependence of Nj(a) on- the basic optical 
properties of the water and its constituents is required in 
order to obtain quantitative information about the _ constituents 
of the ocean from measurements of N, (A) . The remainder of this 
section is devoted to relating N d (Aj to these optical properties. 
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We shall assume that the radiance distribution incident on the 
sea surface is given, and that the radiance is transmitted from 
the ocean to the satellite in a known manner given by a ( a) and 
y (\ ) in equation 2. This reduces the problem to that of taking 
a known downirelling radiance distribution just above the sea 
surface, and computing the distribution (just above the surface) 
of the upwelling radiance. As mentioned above, radiant energy 
interacting with the ocean can be absorbed by water, suspended 
particles, and dissolved organic material commonly called yellow 
substance with absorption coefficients a w , a , a respectively 
and scattered by the water and particles with scattering coef- 
ficients b and (scattering by the yellow substance appears to 
be negligible). The total attenuation, coefficient, c, for these 
interactions is given by 


C = C w + Cp + Cy 


(3) 


where 


C W = a W + t> w> 

c p= a p+bp t (4) 

Cy ‘Qy 


are the beam attenuation coefficients of the water, particles, 
and yellow substance. The scattering is further characterized 
by the phase function, P(0) , which relates to the intensity of 
radiation, dJ (0 ), singly scattered from a small sample volume, dv, -- 
when illuminated by an incidence irradiance, H Q , through 


P(0) = 


dJ(g) 

(H 0 dv)b 


2?r f* P (0) sin0d0 = l . 


(5) 


The total phase function for water and particle scattering is 


✓ 
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'• P ( Q ) = ( b wPw(fl) + bpPp(g)) 

( b w + bp) 

where P w and are the phase functions due to water only 
and particles only respectively. It is convenient to further 
define the forward (F> and backward (B) single scattering 
probabilities by 

Bsl-F 


F=2ir /j" 2 P (8) sind dS 


( 6 ) 


and the single scattering albedo by 


w ° * "c" 


Hence it is clear that 


(7) 


B = 


(bpBp + b w B V v) 
(b p + b w ) 


( 8 ) 


CJ 0 = 


(bw + bp) 


(aw + Qp + ay + bp-f-bw) * 


(9) 


and in general. 


0<<u 0 <| . 


It should be noted that all of the above quantities 
depend on wavelength (A) . The transfer of radiation in the 
ocean is governed by the radiative transfer equation which 
has been discussed in detail by Chandrasekhar (1960) and 
Preisendorfer (1965), Gordon and Brown (1973) have computed 
the diffuse reflectance (upwelling irradiance/incident irradi- 
ance) just above. a flat homogenous ocean as a function of its 
optical properties by a Monte Carlo technique. Using a com- 
bination of the parameters which arise naturally from the 
quasi-single scattering model (Gordon 1973) , Gordon, Brown 
and Jacobs (in preparation) show that the diffuse reflectance 
(R 3 ) can be written 

Rd = O.I79x + 0.05IOx 2 + O.I7IOx 3 Cio) 
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where 


Bm 0 

(l-<u 0 F) 


R d is to first order independent of the distribution of the 
incident irradiance. It should be emphasized that the above 
equation does not include the irradiance specularly reflected 
from the sea surface or white caps, i.e. Rj is the contribution 
to the reflectance from photons which penetrate the surface and 
are multiple scattered back into the atmosphere. Preliminary 
computations (Gordon and Brown unpublished) indicate equation 9 
is also valid for a moderately rough surface. The radiance 
distribution above the sea surface (due to photons scattered 
out of the ocean) is only weakly dependent on , and thus 

the radiance at any viewing angle will in first order vary^with 
o, , and B in the saihe fashion as R d . Hence, it is sufficient 
to study the influence of the optical properties of the ocean 
on R d alone. 

Equation 9 shows that the important oceanic parameters are 
oi 0 and B, and so the observed reflectance spectrum, R^ (A), can be 
explained entirely through a knowledge of «> 0 (a) and B(A). Con- 
versely under optimum conditions, only i» 0 (a) or B(A) can be de- 
duced from R d and then. only. if one of these quantities is al- 
ready known. The situation appears quite depressing when it is 
realised thatw 0 (A) is only imperfectly known even for pure water 
and to our knowledge there are no measurements of B (A) and 
«, (X) for various kinds of suspended particles. Hence at the 
present time it is difficult to interpret «> 0 (A ) and B(a) even if 
both could be extracted from R d (A) measurements. This under- 
scores the necessity of laboratory experimentation to determine 
these optical properties for various ocean constituents such as 
marine phytoplankton, and suspended mineral particles, for ade- 
quate interpretation of oceanic "color M . If the optical prop- 
erties of the constituents are known, then it is theoretically 
possible to determine their concentrations through observations 
Of R d (A) as is discussed in several examples below. 

We shall now examine a particularly simple problem of in- 
terpretation, that of determining the concentration of suspended 
material in the absence and presence of yellow substance. The 
particles are assumed to be nonabsorbing, and their scattering 
coefficient is assumed to be independent of wavelength [cases 
where b p = (const.) A" n with n<l have also been investigated and 
yield results not dramatically different from the n = 0 case). 
These calculations will not apply at all to locations with par- 
ticles containing absorbing pigments such as phytoplankton in the 
water, since Mueller (1973) and Gordon (1974) have shown that 
scattering from such particles varies strongly with wavelength 
near the pigment absorption bands. 
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In order to use equation 10, B (A ) and % (A) in equations 
8 and 9 must be determined. We shall follow Tyler, Smith and 
Wilson (1972) and use' for b w and B w hour g 

computed for sea water continuously filtered for JjJ these 

fPetzold 1972), rather than theoretical computations of these 
Quantities for pure water. We do not feel comfortable with 
this however, on the above basis they have estimated 
abslipUon spectrum of "clear natural water", and since we 
use their a w (x) spectrum this assumption is essential for 
consistency. It is thus assumed that B w - 0.1462 and 
‘ at 550 nm, so 


b w = 0.008m 


b w = 0.008 


' 550 

A 


< 


m 


-i 


To find Br-. Petzold’s measurements of P(0) f ro ™ tongue of 

the Ocean** (T0T0) San Diego Harbor (SDH) j and Off Shore a i 
fornfalSs^ 0 ™ 530 nm ha4 been used Assuming the above for 
R. and b the contribution from water is subtracted yie±uj.n S 
?ge B values listed in Table 1 (« 0 is given in the paren- 
thesil and » 0 for filtered sea water at this wavelength is 

0.136). 


TABLE 1 


Derived Values of B p for the Indicated Locations 


TOTO 

0.0165 (.59) 
0.0158 (.25) 
0.0130 (.26) 


SDH 

0.0186 (.82) 
0.0194 (.83) 
0.0169 (.91) 


OSC 

0.00966 (.59) 
0.00836 (.55) 
0.00836 (.55) 


From Table 1 it is clear that B p not only varies considerably 
from one location to another, but also varies considera y 
for a Tingle region. This , is unfortunate m that for small 
r, i c directly proportional to » B and o 0 ana a are men 
equally important in determining R d for these cases. We 
<;hall assume that Table 1 gives the range of variation of B p 
to be expected in natural waters, however this is probably 
not the case. Furthermore we will assume B p (and B w ) 
independent of wavelength. 

Since a is taken to be zero, and the a w (\) estimate of 
Tyler, Smith p and Wilson is to be used, only a y (A) remains to 
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be considered. This is taken from Jerlov (1968, Page 56) and 
is parameterized by 


a y (X) = a y e ool45 (550-X) 

where \ is in nm. Reasonably high concentrations of yellow- 
substances in the open ocean have ay ~ 5 x 10“^m ' (i.e., 

near Galapagos Islands) however, a v can be much larger in 
coastal regions. y 

Using the above and equation 10, we have computed R(j(X) 
called Reflectance (0,-) in the figures, as a function of 
^p/bw (proportional to the concentration of suspended par- 
ticles) for various values of B p and a . The results are 
given in figures 6 through 10. y ^ 

Figure 6 shows the R^ for \ - 550 nm as a function of 
b p /b w and Bp. The linearity of R^ with b p /b w is to be noted. 
These results clearly demonstrate the importance in knowing 
Bp for a quantitative determination of the particle concen- 
tration from R^ (which is measured at the sea surface') . 

Figure 7 shows spectra of R^j again for various values of 
kp/t>w but now with B p = 0.0165. These spectra have the same 
general shape as those observed by Tyler and Smith (1970) 
just beneath the surface in Crater Lake, but are not in 
quantitative agreement even when the loss due to transmittance 
through the water surface is considered. Basically the 
difference is that quantitative agreement requires too large 
a value of b p /b w . This could be due to the value of B p used, 
the fact B_ is assumed to be independent of wavelength, or 
inaccuracies in their estimated a w (X) on which the present 
calculations are based. In any event the computations can be 
used as a guide for examining ERTS related data. Figure 8 
shows the influence of yellow substance on the reflectance 
for the case with b 
influence of the ye 
region of the reflectance spectrum. All of these results are 
summarized in figures 9 and 10, where the reflectances have 
been integrated over the ERTS MSS channels 4 (figure 9) and 
5 (figure 10) for various B and a y as a function of b r ./b w . 

It is seen that the yellow substance influences the reflec- 
tance in MSS 4 only at very high concentrations (for the open 
ocean) and essentially plays no role in the MSS 5 reflectance. 
Therefore, for suspended particle concentrations , it appears 
that MSS 5 is best with MSS 4 only slightly degraded by the 
dissolved organic materials. It should be emphasized again 
that the above discussion refers only to the case of no 
phytoplankton (chlorophyll) ! 


p/b w - 128 and B p = 0.0165. The major 
llow substance is to depress the blue 
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We now turn to the problem of estimating b p /b w or some 
related quantity from satellite observations. Since the 
above calculations are only for the case of negligible 
phytoplankton (chlorophyll) concentrations we must find 
criteria from which to choose which ERTS channels will satisfy 
this constraint. Also unless only MSS 5 is used we expect 
yellow substances to be important especially in coastal re- 
gions where river runoff, etc., may be considerable. Further- 
more there is the additional problem that the reflectance 
depends on B p which is also unknown, as will as b p /b w . We 
have developed a method which partially overcomes some of 
these problems. Considering the radiances observed in MSS 4 
and 5, we have approximately 


n 4 ■ a 4 N d 4 + y4^SS 4 + N S4 
N 5 = a 5 Nd 3 +75 N SS 5 + N Sg 


where and T i (i = 4, 5) are the fractions of Ndi and N ssi 


(measured at the sea surface) that reach the sensor. It is 
assumed that and are constant over an ERTS Frame (if 
they vary in a known way, their influence is easily accounted 
for and will not be discussed further) . Now from the theory 


^d 4 " k 4 





( 12 ) 


where k 4 and k 5 are essentially independent of b p /b w , but 
depend directly on B p . Taking the horizontal gradient, 
of equation 11 using equation 12 we find 


V H Nj = kjajV H 



+ Xi V H N ss | 


* 


(13) 


since Nc: is nearly constant over a frame (VN~. = 0). The last 

3 1 . 1 „ __ • . t -i - . i - . i-.i 


i^ s -, is the horizontal gradient of the 
the 1 sea surface. This is nearly zero 


term in equation' 13, VfpN 
reflected radiance from t 
everywhere except where the sea state changes dramatically 
with horizontal distance (for example in figure 4). Nearly 
everywhere on the frame then, ^H^ssi - 0 t s° 


V h N 4 k 4 a 4 



Again * 4/ "5 is constant (or slowly varying) over the scene; 
so almost everywhere in the frame, variations in v H N 4 /V H Nc 
are the result of variations in £ 4 / 1 : 5 . Now if 
B p is constant, or if ‘ the "mean" B p is constant in each 
wavelength band, this would imply that essentially only the 
concentration of scattering particles varies over the frame, 
i.e. not the nature of particles or their size distribution. 
In this case, k 4 /k 5 would be constant. Hence essentially 
£ 4 /k 5 will vary, if the nature or size distribution of the 
particles varies over the frame (B p changes), or if the 
yellow substance concentration varies considerably over the 
frame, which would force k 4 to vary independently of ks . 

Thus if we find 


v h n 4 „ 

V h N 5 ' Const (15) 

it is reasonable to expect that only the particle concentration 
changes over the frame. In this case 

VhN 4 orV H N 5 ~v H (particle concentration) ( 16 ) 

Note here that these relations should also apply to scenes 
containing phytoplankton if they are the dominant scatterersj 
but for mixtures of phytoplankton and suspended white particles, 
one would expect equation 15 to be violated over a scene if the 
relative concentrations vary drastically. To reiterate, if 15 
holds, probably only the particle concentration varies over the 
frame, and 16 can be used to measure the gradient of its con- 
centration. Since B p is unknown, the actual concentration is in- 
determinate without ground truth to better than a factor of two, 
since this is the assumed uncertainty in B p . Several examples of 
the use of the above ideas for analysis or ERTS data are pre- 
sented below. 


DISCUSSION 

Equation 16 was applied to MSS 4 and MSS 5 data in 
figure 11. This computer enhanced negative image of the 
Cape Hatteras region shows what appears to be large gradients 
in suspended sediment. The Gulf Stream apparently has 
entrained particles from Rayleigh Bay and is carrying them 
out to sea. Along the scanline shown we have computed 

1 * M * a 17 [n 5 ] + h 
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properties of the individual components with wavelength is 
known . 

* 

As an example of this, the spectra given in figure 12 
were integrated over the MSS 4, MSS 5, and MSS 6 filter 
functions. A series of numerical tests were then made of 
ratios, differences, and sums to see if the three water types 
could be distinguished at the sea surface. It was quite easy 
to distinguish on the basis of such calculations between the 
Gulf Stream waters and the coastal waters, and between the 
coastal waters and the plankton bloom; however it was not 
possible to distinguish between the Gulf Stream and the 
plankton bloom. This suggests, as the theory implies, that 
for these data the ratio test (MSS 4/MSS 5) is not likely to 
be successful in specifying the chlorophyll -a concentration. 
Further, since the sea surface component, N ss , (A) 5 spectrally 
alters N (A), numerical tests (ratios, differences etc.) are 
invalid indicators of oceanic properties. 

Probably the most efficient method of determining the 
concentration of the constitutents in the ocean will be to 
compare theoretical and experimental spectra, adjusting the 
constitutent concentrations in the theoretical spectra until 
agreement is found. This of course requires a basic under- 
standing of the optical properties of the constitutents which 
can be derived only from careful-tn situ and laboratory experi- 
ments. It seems that at the present time much energy and 
money is being expended to try and use optical methods to 
locate and study materials with nearly unknown optical pro- 
perties suspended or dissolved in a medium with only poorly 
known optical properties. This must be overcome before 
significant progress can be made. 

A final interpretation of oceanic observations from ERTS 
is given in figure 13 which is a computer enhanced MSS 6 
negative image of southeastern Florida. The dark lineation 
parallelling the coast in the upper portion of the image is a 
zone of high reflection caused by locally increased N ss along 
the edge of the Florida Current. The increase in surface 
reflectance is probably caused by surface wave interaction with 
the cyclonic boundary and is not bottom influence. This further 
explains the edge effect and supports the discussion on the 
local dependence of VpjN on v fiNss on ly given above. Another 
example of the dominance of N ss is the bright slick areas 
(low N) off the Virginia Key sewer treatment plant. This is 
probably caused by the dampening of the glitter causing 
capillary waves in the oil film associated with the organic 
slick. The slick, which has drifted south past the popular 
Key Biscayne beaches, offers an explanation of the narrow 
lineation off shore in the Florida Straits: a passing oil 

tanker heading south which is pumping her bilges would cause 
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where the least squares values for g=0.38 and h=0.00, and the 
linear correlation coefficient, r, is 0.66. The coefficient of 
determination, r^ , which is the ratio of the explained varia- 
tion from the mean (by the least squares line) to the total 
variation, is 44%. Physically this means that this image may 
not be useful for determining particle concentrations because 
B_, <->o (i.e. ay), or bp/b w may be changing. The extent to which 
this holds true in natural waters is unknown; this will be 
extensively tested in the New York Bight area where turbidity 
measurements are being made by NOAA vessels concurrent with 
ERTS transits. 

The data of figures 7 and 8 can be compared to measure- 
ments of upwelling irradiance. Data given in figure 12 were 
observed using a %-meter Ebert spectroradiometer from 3 meters 
above the surface; these observations were made during the 
time frame of figure 5 and represent the water types shown in 
that image. All spectra were carefully selected to represent 
the same downwelling irradiance, sea state, sun angle, cloud 
cover, and absence of bottom influence. Specular reflection 
due to waves were minimized by preselecting ten spectra with 
similar shapes. After digitizing the records, averages and 
standard deviations (<0 were computed at each wavelength; if 
values exceeded the average by 1«\ they were rejected and a 
new mean computed. Absolute values of the spectra are traced 
to NBS through the 2 -meter integrating sphere at NASA's 
Goddard Space Flight Center. 

As the chlorophyll-a concentration increases, the spectral 
peak shifts to longer wavelengths in similarity with the com- 
puted data in figure 8. This is not to say that the chloro- 
phyll-a and the yellow substance produce quantitatively similar 
results; in general they will not. What it does imply is that 
the broad absorption bands in pigments such as chlorophyll and 
yellow substance at shorter wavelengths (<600 nm) produce 
qualitatively similar effects in the spectra. When other fac- 
tors are equal, increasing the concentration of pigments will 
cause a decrease in the radiance in MSS 4 only. In nature 
however, such as in a plankton bloom or in a river plume, 
increased amounts of chlorophyll are normally accompanied by 
increased particle concentrations (organisms which contain 
the chlorophyll) as well as yellow substances (decay products 
of the plankton) , and increased amounts of "white" suspended 
particles in the river plume case. The radiances in MSS 4 
and 5 will in these cases vary in a manner whieh frill depend 
again on the relative concentration of the constituents , and 
since a and b vary almost independently, the reflectance 
signature, which depends on “ 0 (A) and B(A) is not unique and 
can be unraveled only when the variation of the optical 
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a similar feature on the image. Thus ERTS could he useful in 
coastal waters for such illegal acts which affect 
the nearshore water quality. - 


CONCLUSIONS 

; T £ e spectral properties of the oceanic front associated 
with the Loop Current have been studied by ship and satellite 
observations and by radiative transfer theory. It is seen 
that computer enhancement is required to extract useful 
information from the ERTS data for the ocean scene. The 
current boundary can be detected by changes in the surface 
reflectance, N 55 , as well as the diffuse reflectance, Na» from 
below the surface; N^ however is dependent on both B and 
and thus the spectral interpretation of ocean color requires 
surface truth measurements for meaningful results. Particle 

c a T^ ons » which ca -n delineate currents, can be estimated 
m MSS 5 if the ratio c h N 4 / % n 5 is reasonably constant over a 
scene. MSS 4 is strongly influenced by yellow substance, and 
particle estimation based on these data are invalid in many 
coastal zones. Water mass identification using ratios or 
differences of MSS 4, 5, and 6 data have no validity in either 
theory or observation in the Gulf of Mexico. Finally, it must 
be emphasized that the spectrum of upwelling radiance just 
above the surface is a function of both Nd and N,, and that 
N ss frequently dominates. 
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Time series of the Loop Current in the Gulf of Mexico, Auqust 1972- 
Septemher 1973. The pathlines are the location of the 22°C isotherm 

over the S area eve^fi^f CrU ij e synchronized with ERTS passes 
°; er Iui ar f a J .? ve r y 36 da ^ s * The indicator isotherm was located hv 
expendable bathythermographs from a surface vessel; this isotherm is 

Zn“? ly 20 ki1 ° meterS t0 the n ' ght 0f the ^loni'c edge^f aci ng 
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Figure 2 
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Cross section of the Loop Current taken by ship along a suborbital 
track on the day of ERTS transit. Section shows typical increase 
in salinity and temperature (dashed-radiometric, solid-bucket) when 
crossing into the current. Similarly chlorophyll-a and volume 
scattering function at 45° (436nm light) decrease. The temperature- 
depth section shows the relation between the 22°C at 100 meters and 
the surface frontal zone (at ca. 1330 GMT). Horizontal scale is 5 
hours travel time at 9.8 knots or 90 kilometers. 


16 



1600 


iline is MSS 4, next MSS 5, MSS 6, and 
At scan spot number 950 there is an 
’ 2; this marks the cyclonic edge of 




Figure 4. Negative print of computer enhanced (9<DN_<13; n = 1) MSS 5 image of 
the cyclonic boundary of the Gulf Loop Current, Surface vessel track 
confirmed the location of the current to be the darker shade {higher 
radiance) region on the right hand side of the image (ERTS ID 1065- 
15411), Scanline plot in Fig. 4 horizontally passes through the 
middle of the scene. Horizontal distance across the image is 90 
kilometers. 






Figure 5, negative print of computer enhanced (7<DN<15; n = 2) MSS 5 image of Marquesa Key and Key V ! est (EPxTS 
ID 1099-15293). Change in radiance southwest of Marquesa from dark to light.marks the ship-located 
- boundary between the higher intensity Florida Bay water and the lower intensity Gulf Stream,^ Bottom 
depth is in excess of 100 meters and thus does not contribute to the radiance. Horizontal distance 

across the image is 135 kilometers. 
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Figure 6. Computed reflectance in percent at the sea surface as a function of . the ratio of the particle scat- 
tering coefficient to the water (only) scattering coefficient for 550nm, The value of the fraction 
of backscattered light due to particles (B p ) for each curve is given on the right hand side. 

20 


REFLECTANCE (0,-)% 



wavelength for various values of the ratio of the particle scattering 
coefficient to the water scattering coefficient. Note the wavelength 
dependence in these spectra of changing the particle concentration 
only. 
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Figure 8. Computed reflectance spectra, in percent, at the sea surface. The 
ratio of the particle to water scattering coefficients (b p /b w ) and 
the fraction of backscattered light are 128 and 0,0165 respectively. 
Values of the absorption coefficient due to yellow substance (a v ) are 
listed on the left hand side. Note the shift of the spectral peak to 
longer wavelengths with increasing a y. 
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: igure 9. Computed reflectance in percent at the sea s 
the water scattering coefficients integrated 
0.6pm). Values of the absorption coefficien 
right hand side. Note the linearity depends 

( B p). 
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hgure 10, Computed reflectance at the sea surface, in percent, as a function of the particle to water scat- 
tering coefficient ratio, integrated over the response function of the MSS 5 filter (0.6-0. 7um). 
Valuejjof the absorption coefficient due to yellow substance (a v ) are listed along the right hand 

side. Note the linearity is not strongly dependent on a v but only on the fraction of backscattered 
light, EL. J 



Figure 11. Contrast stretched (4<DN<12; n = 2) negative image of the ocean area offshore of Cape Hatteras 
(CRTS ID 1132-15042). Tlie Gulf Stream can be seen as the bright area to the south of the en- 
trained sediment from the coastal estuaries. The least squares fit of eq. (17) was done along 
a scanline north of the Cape and extending from nearshore, through the suspended sediment and 
into the current. Extensions of this plume were observed for 150 kilometers further east on 
other ERTS images. Horizontal distance across the image is 135 kilometers. 
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Observed upwelling spectral irradiance in the Gulf of Mexico, 
November 1972. The three spectra represent typical observations 
during the time series, and show the shift of the dominant wave- 
length to larger values with increased surface chlorophyll -a. 

The volume scattering function, 645 , is for blue (436nm) light, 
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Figure 13. 


Neqative print of computer contrast stretched (7<DN<15; n = 2) F.RTS 
image of South Florida (ERTS ID 1026-15230). The Florida Current 
can be seen as a line of dark lineation parallel to the coast; bot- 
tom is essentially invisible in this MSS 6 scene. A ship can be 
seen by its characteristic V-shaped wake just offshore of Miami 
Beach. Possibly the Virginia Key sewer outfall area can be observed 
by its low reflectance due to an organic slick. Horizontal distance 


across the image is 90 kilometers. 
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